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Abstract 

We theoretically study the influence of weak interactions on the diffusive expansion of a Bose- 
Einstein condensate in a three-dimensional random potential. For this purpose we develop a per- 
turbative approach and calculate analytically the first-order nonlinear correction to the ensemble- 
averaged atomic density as a function of position and time. We find that interactions renormalize 
the typical diffusion coefficient of the condensate. The magnitude of the nonlinear correction is con- 
trolled by a single dimensionless parameter that includes both interaction and disorder strengths. 

PACS numbers: 03.75.Kk, 71.55.Jv, 05.60.Gg 



I. INTRODUCTION 



Ultracold atomic gases have proven to be a powerful tool to revisit many phenomena of 
mesoscopic physics. In particular, the study of Bose-Einstein condensates (BECs) in random 
potentials allows better understanding of the interplay between disorder and interactions []| . 
The main advantage of experiments with BECs is their great versatility. In particular, the 
strength of the interatomic interaction can be tuned almost at will, making use of Feshbach 
resonances. In the same way, the realization of optical speckle patterns allows an unprece- 
dented control of the parameters of disorder. A prominent experimental example of this is 
the observation of first hints suggesting a Bose-glass phase for strongly interacting BECs 
in disordered potentials [2]. In the last few years, extensive experimental and theoretical 
works have also been carried out on the expansion of BECs in one-dimensional (ID) ran- 
dom potentials 3|, |4|, with particular interest in the phenomenon of Anderson localization 
0, 1^, B]. However up to now, very few works have concerned the expansion of BECs in three- 
dimensional (3D) potentials. In three dimensions the situation is very rich since unlike in 
ID, a critical energy separates localized from extended states js]. Recently, a systematic 
study of Anderson localization of BECs based on the self-consistent theory of Vollhardt and 
Wolfle has been proposed 9|. Somewhat earlier, the diffusive expansion of a BEC in a 3D 
random potential, which predominates when the disorder is weak enough, has been studied 
by Shapiro [loj. In these works however, the role of interactions that come into play during 
the expansion of the BEC in the random potential has been ignored. At long times, when 
the condensate is sufficiently dilute, this assumption is known to be true in the absence of 



disorder 



and seems natural for weak disorder. In the latter case it has been validated 



by numerical simulations in ID 5|. However, even in ID, no consensus on the effect of 
interactions on the expansion of BECs exists ja, and, to our knowledge, there is neither 
analytical nor numerical work about this issue in 3D. 

In this paper we propose the first analytical study of the nonlinear corrections due to 
interatomic interactions to the density profile of a BEC expanding in a 3D random potential. 
We assume that the disorder is weak and neglect Anderson localization. Therefore, we are 
interested in the effect of interactions on the diffusive expansion of the condensate. In the 
absence of interactions the density averaged over an ensemble of realizations of the random 
potential is expected to decay as (D^t) -3 ^ 2 at long times, where is the diffusion coefficient 



2 



of a particle of energy /x, the chemical potential of the condensate [10|| . Surprisingly, we find 
that the first-order nonlinear correction follows the same characteristic decay. Interactions 
thus simply renormalize the diffusion coefficient of the condensate D^, at least at the first- 
order of perturbation theory As expected, we find that repulsive interactions reinforce the 
diffusion unlike Anderson localization which inhibits it. Interestingly, our results show that 
nonlinear effects get stronger as disorder increases. 



II. PERTURB ATIVE TREATMENT 



Consider a weakly interacting Bose-Einstein condensate of N ^> 1 atoms of mass m 
expanding in a three-dimensional random potential V(r), for which we assume a white- 



noise gaussian statistics: V(r)V(r') = K^-k jm 2 ib{Y — r'), where i is the mean free path []i 
and the overline denotes averaging over realizations of the random potential. We treat the 
condensate in the framework of the mean field Gross-Pitaevskii equation 



-^ 2 + V(r)+g\^ 



(1) 



where ip(r,t) is the condensate wave function and g = 47rh 2 a/m is the strength of interac- 
tions, with a the scattering length [11]. Note that the following reasoning is valid for both 
repulsive (a > 0) and attractive (a < 0) interactions. By Fourier transforming Eq. (pQ) with 
respect to time we obtain 



— e 



2m 



V 2 + V(r) 



2^ 



de 2 
2^ 



,(r) = 0, (2) 



where ip e (r) = J dtip(r,t) exp(iet/h). To describe the effect of weak interactions (\g\ <^ 1) 
on the expansion of the condensate we make use of perturbation theory and write ip t {v) = 



V4 0) (r)+^(r), where \6ip e \ < \^'\. A precise condition of validity of this approach will be 
given later on (see section HVBI) . In the following, the superscript ^ will refer to quantities 
free of interactions. Inserting this expression into Eq. (j2J) and keeping only the lowest-order 
terms in 8ip e and g we find that obeys the linear Shrodinger equation 



,(0) 



£ V2 + ^ 



(0)/ 



0. 



(3) 
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whereas the first-order term Sijj e satisfies 



- e - *V> + V(r) 



*(r) = 9 / ^<'C-)V.2'-(r)^l +(1 (r). (4) 

, Z7T ./ 



Eqs. ([3]) and (jl]) form the basis for studying the behavior of the time-dependent ensemble- 



averaged atomic density n(r,t) = \ip(r, t)\ 2 . By expanding n and keeping only the lowest 
terms in Sip we readily obtain 

ra(r,t) ~ n (0) (r,t) + An(r,t), (5) 



where n^(r, t) = |^(°)(r,t)| 2 is the atomic density in the absence of interactions and 



An(r, £) = 2Ke(ip^*(v, t)5ip(r, t)) is the first-order correction to n, with ip^(r, t) and 5^(r, t) 
the Fourier transforms of (r) and 5ip e (r) respectively. The two next sections are devoted 
to the separate study of n^ ' and An. 

III. EXPANSION IN THE ABSENCE OF INTERACTIONS 
A. Diffusive expansion 

In this section we remind a few results concerning the averaged density profile (r, t) 
free of interactions in three dimensions. We assume that the condition of weak disorder 
k^t 3> 1 [3| is fulfilled, where is the wavevector at energy \i. As a consequence, the 
condensate expands mostly by diffusion in the random potential 0, [h| ■ This point will be 
discussed more deeply in section IIIIBI In the absence of interactions, the wave function of 
the condensate at time t and position r is given by 

^°\r,t) = [ ^ [ dVa(r,r')0(r')e-^ /? \ (6) 



2tt _ 

where <fi is the initial condensate wave function and G e is the Green's function of Eq. ([3]). 
From here on, we consider large distances r 3> i and long times t 3> (£/t> M ), where = 
hkn/m is the velocity of a particle with kinetic energy \x. The mean free time tjv^ can be 
regarded as the time separating the ballistic regime t < Ijv^ from the multiple scattering 
one t 3> i/Vfj, where atoms are scattered many times on the random potential. Under these 



assumptions we derive from 
found for example in Ref. [lO]): 



q. ([6]) the following expression for (more details can be 



n • 



d 3 k 



*,*)= I J^MWPe k (r,t). (7) 



Here = h 2 k 2 /2m and P € (t, t) = exp(— |r— r'\ 2 / 4D e t) / (4:7iD e t) 3 ^ 2 is the diffusion propagator 
in three dimensions for a particle at energy e. D e = (l/3)£v € is the diffusion coefficient of 



a particle of velocity v e = ^/2e/m 13]. The momentum distribution of the condensate 
|0(k)| 2 oc (1 — k 2 /2k 2 )H(l — k/y/2kfj), with H(x) the Heaviside step function, is obtained 
from the following expansion scenario j^,!^, 10]: the BEC is generally produced in a harmonic 



trap potential of frequency u. The trap is turned off at some time, which causes the rapid 
expansion of the condensate, driven by strong interatomic interactions. At this initial stage 
of expansion we neglect the effect of disorder. After a time of the order of interactions 
become weak and the momentum distribution is given by ]</>(k)| 2 [14|. We choose our initial 
time t = within this latter stage. 

The calculation of the integral in Eq. (J7]) is straightforward and leads to: 

,<°>M)^/(^y, ( 8 ) 

where f(x) can be expressed through special functions and f(x) ~ 0.04 for z < 1. As 
pointed out in 10], the density profile is driven by a single parameter = (l/3)v ll £, which 
can be regarded as the typical diffusion coefficient of the condensate. We plot the density 
profile in Fig. [T]as a function of time for three different distances r. n^(r, t) exhibits the well 
known decay in t -3 / 2 in the long time limit t ^> tamvai, where t arr ivai ~ r 2 /D^ ~ {l/v^){r / ' t) 2 
is the "arrival "time at which the density reaches the maximum. 

B. Anderson localization 

It must be noted that Eq. (jSJ) describes a purely diffusive expansion and does not include 
the effects of Anderson localization. In Ref. the authors found oc 1/t for t > £ arr ivai- 
A question then naturally arises: to which extent the diffusive result (jSJ) is correct? In Ref. 
9j the density profile has been studied for moderate values of k^i, that is in the regime of 
strong disorder. When k^i is increased, the 1/t decay is expected to cross over to the l/t 3/2 
decay typical for diffusion. To illustrate this cross-over we plot n^°'(r, t) including localization 
effects in Fig. [2] for three values of k^i and for r = 10£ (three upper curves, orange, green 
and blue online). The lower curve in Fig. [2] (red online) is obtained from Eq. (jSJ), i.e. by 
neglecting Anderson localization effects. For k^i = 3, a significant fraction of atoms of the 
condensate are localized. At very short times, all diffusive atoms have already flown away and 
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FIG. 1: (color online). Ensemble-averaged atomic density n^°'(r,t) of a BEC expanding without 
interactions in a 3D random potential, in the regime of weak disorder k^l 3> 1, as a function of 
time and for three different distances r = |r| from the initial location of the condensate. The 
dashed line is a 1/t 3 ^ 2 asymptote. The time scale is in units of ^arrival such that all maxima fall in 
the same abscissa. 

n^ -* oc 1/t, in agreement with Ref. For weak disorder {k^i = 30inFig$2§ } the situation 
is very different: most atoms of the condensate propagate diffusively. Therefore the atomic 
density exhibits a 1/t 3 / 2 decay characteristic for diffusion. Whatever k^t however, at very 
long times the atomic density saturates at a constant value n[^ c ~ h(k fl £)(N/r 3 )(£/r) 1 ^ , 
where v is the critical exponent of the localization transition, and hik^l) oc 1/(&; M £) 3 for 
k^l 3> 1 and h{k^t) oc const for k^l « 1 |. We can estimate the time after which the 
diffusive approach breaks down by requiring that the saturation value n^ c is smaller than 
the limit nJJ ~ N/^D^tf/ 2 of Eq. © at long times t > t arr ivai- This yields 

tioc~(^)(r/£) 2{1+1/3 ^(M) 2 , (9) 

as the maximum time at which diffusion model can be applied. When k^i 3> 1 and r ^> i, 
t\ oc is always larger than t a mvai and the diffusion model is hence valid in a broad time 
interval. To sum up, the condensate behaves diffusively (n^ oc i~ 3 / 2 ) until t = t\ oc , and 
beyond localization effects start to "freeze" the density profile at = n[^ c - The conditions 
k^£ ^> 1 and t <C t\ oc will be assumed throughout the remainder of the paper. 
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FIG. 2: (color online). Ensemble-averaged atomic density n^°\r,t) of a BEC expanding without 
interactions in a 3D random potential, for three different values of the disorder parameter k^l and 
for r = 10^ (orange, green and blue upper curves). Here unlike Fig. [TJ the effects of Anderson 
localization are fully taken into account. For moderate values of k, t £, the curves exhibit first a lit 

n 

decay and then saturation at some constant value, in agreement with Ref.[9]. When k^i ^> 1, the 
atomic density exhibits a l/i 3//2 decay, before saturating at very long times as well. The red lower 
curve is obtained from Eq. ([8]), i.e. by neglecting Anderson localization effects. 

IV. NONLINEAR CORRECTION TO THE ATOMIC DENSITY 

In this section we are interested in the first-order nonlinear correction An to the density 
profile [see Eq. (jSJ)]. We first discuss the effect of interactions on the mean free path, and 
then develop a diagrammatic approach for calculating An. 



A. Mean free path 

It should be noted that the mean free path i is in principle modified in the presence of 
interactions. However in the following we will neglect the nonlinear corrections to I. This 



approximation is legitimate in 



;he limit of long times, as we now show by making use of a 



criterion initially developed in 15] for monochromatic classical waves. 



The mean free path is weakly affected by interactions provided that the scattering on the 
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effective potential g\i/)(r, t)\ 2 is weak as compared to the scattering on the random potential. 
In other words, the mean free path £^ associated with the interaction potential should be 
much larger than the mean free path £^ = £ associated with the random potential. i^ 1 ' can 



be evaluated from the short-range "Ci" correlation function of the density 5n(r,t)6n(r',t) 
which has been recently studied by Henseler and Shapiro 16]. We neglect here the con- 
tribution due to the long-range "C 2 " part of the density correlation function. This would 
be questionable for a condensate propagating in a quasi-one dimensional waveguide where 



C*2 was predicted to dominate at long times 171 ] ■ However the situation is different in an 
unbounded 3D random potential where we expect C2 to decay with time. Indeed, let us 
consider a quasi-monochromatic wave pulse emitted from some point in an unbounded three 



dimensional disordered medium. A straightforward calculation similar to that of Ref. [17 ] 
yields C2 oc (l/t 2 )/(k£) 2 in the limit of long times. C2 is therefore negligible with respect to 

c " n 

From 16[] we write 5n(r,t)5n(r',t) ~ (2Tr£/k 2 )n^°\r, t) 2 5(r — r'). For times t ^> t a mvai, 
this leads to ~ (h 4 Tr/m 2 )( y k 2 /27r£)(D^t) 3 /( y g 2 N 2 ). The condition > £ then gives 



t > ti = p 



2/3 



(M) 1/3: 



(10) 



where we introduced a dimensionless parameter p = aN/£^k^£ that includes both inter- 
action and disorder strengths. At this point, it is sufficient to keep in mind that in the 
weak disorder limit and under conditions of typical experiments, p< 1. For example with 



the data from Refs. 



and 



6j we find p ~ 10 2 . Therefore, the characteristic time t\ is 
much smaller than £jv^ such that the criterion ( TlOl is automatically satisfied for long times 

t > w 



B. Diffusion coefficient 



We now consider An(r, t). From Eq. (TJJ we obtain 

<ty(r,t) = J^J dVG f (r,r')%')e^ £l/fi , 



where S(r') is the right-hand side of Eq. 
we obtain 



evaluated at r'. Combining Eqs. 



(11) 
and (fill) 
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An(r,t) = 2Re V (0) *0, t)5ip(r, t) 



2Re 



! Yldejdhjdh'e^^-^Kir, t; r', {r,}, {e^}) 

^ 3=1 



x0( ri )0*(r 2 )0(r 3 )(/»*(r 4 ) 



(12) 



where the six-point kernel is given by the connected part of a product of five Green's 
functions, averaged over disorder: 



K = G ei (r,r')G* 2 (r,r 2 )G e3 (r',r 3 )G: 4 (r',r 4 )G ei . 



(13) 



The problem therefore reduces to finding the largest contribution to the kernel K in the 
weak disorder limit k^i ^> 1. This contribution is given by the diagram depicted in Fig. 
|2k>. The diagram describes interaction of two matter waves at a point r'. Diagrams of this 



type have been introduced some time ago 



of nonlinear optics of disordered media 



'or monochromatic classical waves in the context 



181 ] . In our case the situation is more general since 



K contains a product of Green's functions at five different energies. Writing down explicitly 



a) 



?3 



T4 



^3 



ri- 



^1 — + 



b) 



1*2 



€2 

FIG. 3: a) Diagram for the kernel K giving the first nonlinear correction to the atomic density. 
The solid and dashed lines represent Green's functions G ti and G* , respectively, where ej = e 4 , e 3 
or e\ — e 3 + e 4 and €j = e 2 or e 4 . The two parallel G lines connected by dotted "ladders" symbolize 



averages of products of two Green's functions, G ei G* . b) Diagram J of Eq. (|14p . 
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the diagram of Fig. [3^, we obtain 



An(r,f) = 2Re 



(-i/R)(ei-e a )t, 



(r 1 )^(r 2 )0(r 3 )0*(r 4 ) 



xG £l (r,r')G* (r,r') x G £3 (r', r 3 )G* (r' ; r 4 ) x G ei _ £3+e4 (r', n)G* (r', r 2 ) 



,(14) 



where p = mk^/^^h 2 is the density of states at energy /x and J = J d 3 ~k/ (27r) 3 G(k)G*(k) = 
—i(£/v fl )(Air 3 p 2 £/hk 2 l ) is the diagram depicted in Fig. [3b. Now we introduce new variables 
SI = €i — e 2 , f2' = e 3 — e 4 , e = (ei + e 2 )/2, e' = (e 3 + e 4 )/2, and make use of results of Refs. 

urn 



d 3 r 3 d 3 r 4 G* (r 3 ,r')G* (r 4 , r')0(r 3 )0*(r 4 ) = -- 



d 3 k' 



P £ ,(r',^)lm^(k')|0(k')| 2 , (15) 



n J (2tt) 3 ' 

where P e (r,Sl) = exp(— |r| a/ —iSl/HD e ) / (4wD e \r\) is the Fourier transform of the diffusion 
propagator at energy e. A similar relation holds for the integrals over ri and r 2 . In the weak 
disorder limit k^£ ^> 1, ImG e '(k') can be approximated by — 7r5(e' — ey) flHQ. ths 
finally leads to 

An(r,t) = 2Re 



gj vf 



d 3 v' 



d 3 k f d 3 k' 



dSl / dtt'e- int/H 



pn) 3 ph£ 2 J J (2tt) 3 J (2tt) 3 
k)H0(k')| 2 Pe fc (r',fi - 0')^( r, » fi/ )^+nv 2 (r - r', SI) 



(16) 



The main difficulty of Eq. (TT6"j) lies in the four coupled integrals over energies and mo- 
menta, which make the exact calculation of An complicated. However, after some algebra 
and in the limit of long times, all integrals can be performed. The main lines of this calcula- 
tion are reported in Appendix IA1 We find that for t 3> t 2 = (£/v /i )(r/£) 4 (fc M £), An becomes 
position independent as n^ , and 

N 



An(r,£) 



-P 



{D,tf/ 2 - 



(17) 



It is quite remarkable that the first-order correction An has exactly the same time depen- 
dence as n^ ) at long times. Comparison of Eq. (fT7j) to 



(0) 



N/(D fl t) 3/2 provides the 



criterion p < 1 as a condition of validity of our perturbative approach. As discussed in 



section IIV A[ this criterion is usually satisfied in typical experiments, which validates the 
assumption of negligible interactions during expansion made in previous works [lol. 

Few comments are in order. First, the negative sign in Eq. (fT7j) indicates that repulsive 
interactions (p > 0) tend to reinforce the diffusion process and the contrary for attractive 
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Diffusive regime 



Region of validity [ Localized 
of Eqs. (17) and (18) regime 



regime 





time 



v 



FIG. 4: Ordering of the different time scales that we use throughout the paper, t/va is the mean 
free time, t\ is the time after which the mean free path is not affected by interactions, tarrival is the 
time at which the atomic density free of interactions is maximum, ti is the limit of long times for 
nonlinear corrections to the atomic density, and t\ oc is the time after which Anderson localization 
cannot be neglected anymore. In this scheme, the limits of weak disorder k^l S> 1 and large 
distances r»f are assumed. 

interactions, which might have been expected. This appears more clearly if one writes the 
total density profile as 



where we have defined an effective diffusion coefficient D e R = D^/{1 — p) 2 ^ 3 taking interac- 
tions into account. For repulsive interactions for example, p > and therefore D c g > D^. 
Thus repulsive interactions seem to compete with localization, even though calculations of 
higher-order diagrams or a complete study of interactions within the localized regime would 
be necessary to confirm this effect. For extremely weak disorder (k^£ — > oo) p vanishes and 
D e ff ~ whatever the interaction strength. 

Another important remark is that the magnitude of the correction term (fT7|) can be 
modified by increasing the strength of disorder independently of the strength of interactions, 
or vice versa. In particular, if we extrapolate our calculation to k^i ~ 1, we find that even 
weak interactions may significantly affect the expansion of the condensate. In this case, the 
scattering length a must be much less than £/N to justify the neglect of interactions. This 
refines the condition of validity for the results of Ref . |9|] , where interactions were assumed 
weak, but no precise condition of weakness was given. 

Our result (|T7|) is only valid in the time interval [^^loc]- F° r weak disorder, t 2 <C ti oc , 
which gives a broad range of validity to Eq. ffTTl) . For the sake of clarity we summarize all 
the time scales that we introduced above in Fig. HI 




(18) 
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V. CONCLUSION 



In conclusion, we developed a perturbative technique to study interaction corrections to 
the density of a Bose-Einstein condensate expanding in a three-dimensional random poten- 
tial. Our results apply to the case of weak disorder. We have shown that at the first-order 
of perturbation theory and at long times, interactions simply renormalize the effective dif- 
fusion coefficient of the condensate. Repulsive interactions between atoms accelerate the 
diffusion process whereas attractive interactions slow it down. In addition, we have found 
that the effect of interactions on the average density is controlled by a single parameter 
p, depending on both disorder and interaction strengths. Remarkably, the corrections to 
the density profile tend to increase when disorder becomes stronger, at a fixed strength of 
interactions. Our main result (IT7|) has been derived under the assumption of a white-noise, 
uncorrelated random potential. The generalization to correlated potentials used in experi- 
ments is straightforward. For example for a speckle potential of amplitude Vr and standard 
deviation <xr, one simply has to replace the mean free path by k^h 4 /7cm 2 V^<j R [h]]. The 
calculation of higher-order interaction corrections would be of considerable interest in order 
to validate our main conclusions, even if this rises important yet purely technical difficulties. 
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APPENDIX A 

In this appendix we give the main lines of the derivation of Eq. (TiTj) from Eq. ffTBT) . The 
first step consists of transforming both integrals over k and k' into integrals over energies e 
and e' by making use of e = h 2 k 2 /2m and e' = h 2 k' 2 /2m. Expressing explicitly the diffusion 
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propagators, Eq. (fISl) can be written as 



An(r,t) = 2Re 



g J v; fm y [°° r ,2 dr , i sin( )d / de 



(2ixf phi 2 \h?J J 

xv / ii 7 |0(e)| 2 |0(e , )| 2 exp 
exp(— iQt/h) 



•x p2e p2e' 

de' I dtt i dVL' 



r — r 



hD, 



e+Cl'/2 



-in' /-i(n-fi') 
tld7 + 



Here the bounds of integration over Q and f2' are obtained by noting that e ± f2/2 and 
e' ± fi'/ 2 are originally energies (see the main text) and therefore positive quantities. As 
a consequence, \Q\ < 2e and \Q'\ < 2e'. The energy distributions |0(e)| 2 and |0(e')| 2 are 
deduced from momenta distributions |0(k)| 2 and |</>(k')| 2 (we keep the same notation <ft for 
clarity). We have for example: 



(Al) 



■ 2 ( 

1 (v^M 3 2 ^ 



' "\fffl- ' 



(A2) 



with H the Heaviside step function. In Eq. (1A2|) the prefactor has been obtained by requiring 
the conservation of the total number of atoms N = J c/ 3 k|0(k)| 2 /(27r) 3 . The integral over 9 
is readily performed: 



sin(0)d0- 



r — r 



xrr 



e -x\r-r'\ _ e ~x(r+r') 



(A3) 



where x 



hD, 



e+VL'/2 



-. Eq. (IAip then reduces to: 



An(r,t) = 2Re 



gj w 2 (my 



_(2tt) 9 P W? \h?J 




rD t+n i /2 D e iD t 

where we introduced a lower cutoff A ~ I in the integral over r' . This cutoff is needed 
because of the breakdown of diffusion theory at small length scales. The next step consists 
of calculating the integral over r' . Before doing so it is convenient to introduce the new 
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dimensionless variables u = Qt/2h, v = Q't/2h, p = et/h, q = e't/h, r = k^£t/(£/v^), 
p = (r/£)^3h~i/V^, p> = (r/^^V/v^/r 1 / 4 and A' = (A/£) ^3 k/fV^/r 1 ^. Using 
D e = (£/3)^/2e/m we obtain 



An(r,t) = Re 



-iC 



00 dp' r 

A' P' Jo 



V 



dp j d(/ j du I dv-^_ 7 ,., 



1 T l ' A 



VP + v ( x 
-2iu \ 



exp — p 



-2iv 

7^ 



-2i(u — v ) 



VP 



x ^exp 
exp(— 2iu) 

Vp + v 



-1/4 



-2iu 



VP + v 



exp 



-1/4 



-2iu 



VP + v 



(AS) 



We can now perform the 



where the prefactor C ~ (gmN 2 k ll €) / '(hi) ~ (aN 2 k i 
integral over p'. In the limit p <C r 1 / 4 which corresponds to long times t t 2 (t 2 is defined 
in the main text), this integral is 



df/_ 

A' p' 



exp 



-2iu 



VP + v 



exp 



-2iu 



VP + v 



x exp — p 



-2iv 



+ 



-2i{u — v) 



T 



1/4 



-A' 



-2iw 



, P 



—2iu 



(A6) 



where the last equality results from p/r 1 / 4 ^> A'. Inserting Eq. (1A6|) into Eq. (1A5[) we 



obtain 



An(r,t) = 2Re 



C 



T 



7/2 



dp / dg / du dv 



exp (—2m) 
y/p + V 



The four remaining integrals can be readily performed. We finally obtain: 

aN 2 1 



MM)- -572 



(A7) 



(A8) 



which is Eq. f lT7|) of the main text. 
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